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ABSTRACT 

Stochastic acceleration of charged particles interacting resonantly with a tur- 
bulent magnetic field in a small accretion torus appears to be the likely mecha- 
nism responsible for much of Sagittarius A*'s mm and shorter wavelength spec- 
trum. The longer wavelength radiation is produced at larger radii by electrons 
either diffusing from smaller scales or accelerated in situ. An important pre- 
diction of this model is the ejection of a significant fiux of relativistic protons 
from a magnetic-field-dominated acceleration site into the wind-shocked medium 
surrounding the black hole. Recently, several air Cerenkov telescopes, notably 
HESS, have detected TeV emission from within V of the Galactic Center, with 
characteristics hinting at a pp-induced pion decay process for the 7-ray emis- 
sion. Given (i) that we now know the size of this acceleration region, where 
Sagittarius A*'s 7-mm wavelength emission originates, and (ii) that we can now 
map the wind-injected ISM within ~ 3 pc of the nucleus using the diffuse X-rays 
detected with Chandra, it is feasible to test the idea that protons accelerated 
within ~ 20 Schwarzschild radii of the black hole produce the TeV emission far- 
ther out. We show that the diffusion length of these particles away from their 
source guarantees a majority of TeV protons scattering about once within ~ 3 pc 
of Sagittarius A*, and we demonstrate that the proton power (~ 10^"^ ergs s~^) 
produced in concert with the 7-mm radio emission matches the TeV luminosity 
well. The particle cascade generated by the pp scatterings produces GeV 7-rays 
from TT^ decays, and bremsstrahlung, inverse Compton, and synchrotron emission 
at longer wavelengths from secondary particles. We compare these with current 
measurements and demonstrate that GLAST will detect this source during its 
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one-year all-sky survey. This model explains why the TeV source is unresolved, 
yet does not vary on a time scale of a year or less, and it also accounts for the 
high-energy emission while retaining consistency with Sgr A*'s well-studied cm 
and mm characteristics. 

Subject headings: acceleration of particles — black hole physics — Galaxy: center 
— plasmas — turbulence 

1. Introduction 

Sagittarius A*, the compact radio source associated with the supermassive black hole 
(of mass M ^3Ax 10^ M©; Schodel et al. 2003, Ghez et al. 2003) at the Galactic Center, 
has a stratified emission structure in which the higher frequency emission is produced from 
progressively smaller volumes (see, e.g., Zhao, Bower, & Goss 2001; Genzel et al. 2003; 
Baganoff et al. 2001; Meha, Jokipii, & Narayanan 1992). The variable component of the 
near-IR (NIR) to X-ray emission fluctuates on a time scale of a few minutes (Porquet et al. 
2003; Ghez et al. 2004) and exhibits a quasi-periodic modulation during some long duration 
flares (Genzel et al. 2003; Belanger et al. 2005), hinting at an origin within an accretion 
torus with size equal to a few Schwarzschild radii (r^ = 2GM/c^ ~ 10"*^^ cm; Liu, Melia, 
& Pctrosian 2006). The spectral (Falcke et al. 1998) and polarization (Aitken et al. 2000; 
Bower et al. 2005) characteristics of the mm/sub-mm radiation suggest that these too are 
produced within such a region (Meha, Liu, & Coker 2000). The cm and longer wavelength 
emission, on the other hand, has to be produced on progressively larger scales to avoid being 
self-absorbed (Liu & Melia 2001). Recent observations by Bower et al. (2004) show that 
indeed Sagittarius A*'s photosphere has a radius ~ 24r5 at 7 mm. 

In earher work (Liu, Petrosian, & Melia 2004, hereafter LPM04), we showed that the 
stochastic acceleration (SA) of electrons interacting resonantly with plasma waves or turbu- 
lence generated via an MHD dissipation process (see, e.g.. Miller, LaRosa, & Moore 1996; 
Hamilton & Petrosian 1992; Petrosian & Liu 2004), can produce the nonthermal particles 
emitting much of Sagittarius A*'s spectrum (see, e.g., Zhao, Bower, & Goss 2001; Genzel et 
al. 2003; Baganoff ct al. 2001). This mechanism, originally invoked to explain how magnetic 
reconncction can energize particles during solar flares, can also accelerate electrons in the 
magnetized Keplerian accretion torus within ~ 5 — 10 of the black hole via a turbulence 
induced by the magneto-rotational instabihty (Balbus & Hawley 1991; Meha, Liu, & Coker 
2001; Hawley & Balbus 2002). The radiation emitted at the acceleration site can explain 
the mm/sub-mm and shorter wavelength observations of this source in both the quiescent 
and flaring states. 
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An interesting prediction of this model is that, under most circumstances (see Fig. 1 
in LPM04), the time required for energetic electrons to escape from the acceleration site 
is comparable to the acceleration and synchrotron cooling time scales so that there would 
be a significant outfiux of accelerated particles. The nonthermal "halo" or jet formed by 
the escaping electrons may contribute to the cm- wavelength emission farther out (see, e.g., 
Liu & Melia 2001; LPM04; Zhao et al. 2004). In concert with the acceleration of electrons 
within ~ 5 — lOrs of the black hole, protons are also energized. But as we shall show in this 
paper, their acceleration efficiency in this "gas-dominated" region (with Alfven velocity va 
much less than the speed of light c) is not sufficient to produce an observationally important 
flux away from Sagittarius A* itself. 

However, particles continue to be accelerated outside of the inner torus, as evidenced, 
e.g., by the intensity and scale of the 7-mm emission region. We will therefore here also 
consider the in situ SA of particles in a larger "magnetically" dominated {va — c) volume 
(extending out as far as ~ 20 — 30r5), and demonstrate that the conditions required to 
account for Sagittarius A*'s 7-mm emission are in fact equally suited to the rapid acceleration 
of protons that escape and interact with the ambient medium surrounding Sagittarius A* to 
produce 7-rays. 

Recently, the Galactic center (GC) has been identified by three air Cerenkov telescopes 
at ~ TcV energies: Whipple (Kosack et al. 2004), CANGAROO (Tsuchiya et al. 2004) and, 
most recently and significantly, HESS (Aharonian et al. 2004). The Hcgra ACT instrument 
has also put a (weak) upper limit on GC emission at 4.5 TeV (Aharonian et al. 2002) and 
the Milagro water Cerenkov extensive air-shower array has released a preliminary finding 
of a detection at similar energies from the region defined as I e {20°, 100°} and |6| < 5° 
(Fleysher 2002). We examine here the likehhood of these TeV photons being produced by 
the relativistic protons escaping the SA site near Sagittarius A*, and concentrate on the 
HESS results, though it should be noted that all of these observations lend crucial support 
to the notion that acceleration of particles to very high energies is taking place near (or at) 
the black hole. Observations pertinent to these high energy processes are described in § 2. In 
§ 3 we discuss the acceleration of electrons and protons by turbulence and the observational 
constraints on the magnetic field and gas density of an acceleration region roughly 20r5 
in size. The process of production of TeV photons is described in § 4 and a summary is 
presented in § 5. 
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2. Contextual Background 

HESS detected a signal from the Galactic center in observations conducted over two 
epochs (June/July 2003 and July/August 2003) with a ---^ 6 — 9a excess evident over this 
period. These data can be fitted by a power law with a spectral index 2.21 ± 0.09, with a 
total flux above the instrument's 165 GeV threshold of (1.82 ± 0.22) x 10"^ m'^ s"^ (there 
is also a 15-20% error from the energy resolution uncertainty). 

Intriguingly, the HESS data are difficult to reconcile with the EGRET GC source 
3EG J1746-2851 (generally ascribed to pp-induced neutral pion decay at GeV energies; 
Fatuzzo & Meha 2003). In a re-analysis of select data from the 3EG catalog (Hartmann 
et al. 1999), Hooper & Dingus (2002) have shown that the GC is excluded at the 99.9% 
confidence hmit as the true position of 3EG J1746-2851. The HESS source, on the other 
hand, is coincident within ~ 1' of Sagittarius A*, though its centroid is displaced roughly 10" 
(corresponding to about 0.4 pc at a distance D = 8 kpc) to the East of the GC. In addition, 
an extrapolation of the EGRET spectrum into HESS's range over-predicts (by a factor ~ 20) 
the TeV 7-ray flux of the GC source. Recently, Crocker et al. (2005), and Aharonian & 
Neronov (2005a, 2005b), among others, have undertaken a detailed examination of the pos- 
sible theoretical consequences of these new high-energy observations of the Galactic center, 
and inferred that we are apparently dealing with two separate sources. Perhaps both are 
associated with the SNR Sagittarius A East, but given its angular proximity to Sagittarius 
A*, the TeV source could be associated with the black hole itself. 

What is not yet clear, however, is how the TeV 7-rays are produced. Models that invoke 
pp-induced neutral pion decays close to the black hole must contend with the collateral 
emission at other wavelengths (Markoff, Melia, & Sarcevic 1999; Aharonian & Neronov 
2005a) due to the secondary particles produced in the ensuing cascade. And because the 
inner ~ 10 r^-emitting region of Sagittarius A* has received extensive scrutiny over the years 
(see, e.g., Melia, 1992; Falcke et al. 1998; Bromley et al. 2001; Yuan, Quataert, & Narayan 
2003; Melia & Falcke 2001; Melia 2006), non-hadronic models face similar difficulties in 
satisfying the constraints imposed across Sagittarius A*'s broadband spectrum — from radio 
waves to X-rays and 7-rays (but see Atoyan & Dermer 2004 for a description of the "plerion" 
scenario, which avoids this complication). In the model we are proposing here, however, pp 
scattering events are quite rare in the acceleration site since the ambient proton density is 
very low. Instead, most of the accelerated protons diffuse out of the system and scatter with 
the ISM surrounding the black hole. 

Our understanding of the circum-black-hole environment at the GC has improved con- 
siderably in recent years, thanks to observations with Chandra and detailed hydrodynamic 
modehng of the diffuse X-ray emission it detected within r ~ 10" of Sagittarius A*. The hot 
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plasma producing this radiation, estimated at 7.6 x 10^^ ergs s ^ arcsec ^ in the 2 — 10 keV 
band, has an RMS electron density {n) ~ 26 77^^^^ cm~^ and a temperature k^T f^i 1.3 keV, 
with a total inferred mass of ~ 0.1 MQirj^J"^ (Baganoff et al. 2003), where r// is the volume 
filling factor. We now understand that this hot gas is produced via collisions of winds of stars 
in the central cluster injected into the mini-cavity with a total mass rate of M ~ 3 x lO'^M© 
yr~^ (Rockefeller et al. 2004). In the most recent census, some 25 bright, young wind sources 
have been identified and their mutual interactions account for the entire X-ray flux observed 
with Chandra. For a typical wind velocity v — 1000 km s"-*^, the average (pre-shock) gas 
density in the central ~ 1 pc region is (no) ~ 200 cm~^. We therefore have a good handle on 
the physical conditions in this region. In addition, this medium appears to be threaded by 
a rather strong magnetic field ~ 0.1 — 1 mG (Yusef-Zadeh et al. 1996), through which the 
escaping relativistic protons must diffuse. Note that this measured magnetic field is roughly 
in equipartition with the gas. 



3. The Acceleration Site 



To date, the most rehable measurement of Sagittarius A*'s intrinsic size has been made 
with closure amplitude imaging using the VLBA (Bower et al. 2004). At a wavelength 
of 7 mm, the source is confined to the inner ~ 24 Schwarzschild radii, though the size is 
variable on a time scale of hours to weeks and sometimes reaches ~ 60 r^. Correspondingly, 
Sagittarius A*'s average flux density at this wavelength is 1.5 Jy (e.g., Falcke et al. 1998). 
These set strict constraints on the SA model, because the acceleration must not only be 
confined to a region of known size (~ 24 r^), but must also produce a particle distribution 
that accounts for the measured radiative output at this wavelength. 

A complete treatment of SA by waves requires a solution of the coupled kinetic equations 
for the pitch-angle-averaged particle distribution N{E,t) and the wave spectrum >V(k): 
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where = 7— 1, 7 and k are, respectively, the particle kinetic energy in units of the rest mass 
energy, Lorentz factor and the wave vector. The loss rate, the source and escape terms for the 
particles are given by El, Q and N{E)/Tcsc, respectively. The terms on the right-hand-side 
of Equation (2) represent the wave generation, damping, cascade, and leakage processes. The 
energy diffusion coefficient Dee{E) and the function A{E) = {Dee/ E){2 — 7~^)/(l + 7^^) 
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depend on W(k), and the wave damping rate r(k) depends on N{E). That is how the 
equations are coupled. 

SA of particles is dominated by the transit-time damping and cyclotron-resonance pro- 
cesses. Recent progress in studying the scattering and acceleration of charged particles by 
MHD turbulence indicates that fast mode waves play the dominant role (Yan & Lazarian 
2002). Fast wave turbulence is isotropic and has a power-law spectrum with an index of 
g = 1.5 in the inertial range (Cho, Lazarian & Vishniac 2002). At small scales, the waves 
are subjected to collisionless damping, which depends on the wave propagation angle with 
respect to the magnetic field 6 (Stepanov 1958; Ginzburg 1961; Petrosian, Yan, & Lazarian 
2005): 

r(k) c ^7r/3/)V2^^^^ , for 1»P,»S, (3) 
z COS a 

where u, S = me/nip. and j3p = AnPgas/B^ are the ratio of the proton to electron en- 
ergy density, the electron to proton mass, and the gas to magnetic field pressure, respec- 
tively. The wave spectrum will cut off when this damping rate dominates the cascade rate 
'''cas(^) — '^A/turb(fcmm^/27r)^/^, whcrc 2n / k^^i^ gives the turbulence injection length scale and 
/turb = [Stt J W{k)dk]/ B'^ < 1 is the ratio of the turbulence to magnetic field energy density. 
Most of the waves at small scales therefore propagate along the magnetic field fine with the 
propagation angle 9 < 9ci{k) — {2f'^^^^k^i^/ t^'^u^ PpSky/^ . To study particle acceleration by 
such turbulence, one may represent the latter as parallel propagating waves with spectrum 
W oc k-^/'^el oc A;-2. 

Almost all of the observed emission from the direction of Sagittarius A* is produced by 
nonthermal relativistic particles (£■ ~ 7 » 1 and A ~ 2Dee/E). As such, the characteristic 
interaction time Tp, the pitch-angle-averaged scattering (tsc = ([1 — /i^]^ / 0^^)) and accelera- 
tion (Tac = E'^/Di^e) times, and the escape time, can all be expressed in units of rtr = R/c, 
the light transit time for an acceleration region of size R (Dung & Petrosian 1994; Petrosian 
& Liu 2004): 

^ ^ 1 f ^^^A 3 f ^ (4^) 

" 7r(g - l)/turb^^ V ^ ) ' - l)/turbi?fcmin V ^ J ' 

Ttr f -jck, ^ 



T, 



/Si ' /Si - l)/turbi?A;^in V ^ 

n-. 



2-q 



2{(3l + l)Tp_,_^ ^ 2.6(/?i + 1) Ttr f ick,,,„, ^ 



' sc 



o.8(g - i)/turbi?^„.i. ( ^%^y '+1.4 



Ttr, (7) 



n 

where /i, -D^//, and Q are, respectively, the cosine of the particle pitch angle with respect 
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to the large scale magnetic field B, the pitch angle scattering rate, and the nonrelativistic 
gyrofrequency of the particles, and 27r/A;inin < R. The (dimensionless) Alfven velocity is 
defined as 

~ (47rnmpc2)V2 ^-"^ \iG ) llcm-s; " [1 - (^;a/c)2]V2 " 

Note that, for /9a > 1, the phase velocity of the waves is va/c = w/{ck) ~ /9a/(1 + 0aY^'^- 

For q = 2, the acceleration and escape time scales become energy independent. When 
the loss processes are unimportant (i.e.. El <ti A), Equation (1) can be reduced to 
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whose steady-state solution is N ex E p above the particle injection energy, where the 
spectral index p is determined by the ratio of Tac to Tesc' 

9 r 

•-' ' ac 



P-[l + ^J -0.5 > 1.0. (10) 

There are four primary parameters in the SA model (LPM04): R, /turb-R^min, n and 
B. For Sagittarius A*, i? ~ 20r5 at the wavelength of interest (i.e., 7 mm). The observed 
7-ray spectral index and the radio to IR spectrum set strict constraints on the other three. 
To produce a power-law spectrum of relativistic protons with p = 2.2, which gives rise to a 
7-ray spectral index of the same value via pp scattering processes. Equations (10) and (7) 
show that Tac — 5.0 Tesc = (4/turb-R^mm + 7)Ttr. Combining this with Equation (6), one has 

2.6 ^ 



U[/turbi?A;^in]' + 7/turbi?A;^in-2.6j ■ ^^^^ 

Since /turb < 1, this sets a constraint on /3a- 

Next we consider the constraints imposed on B and n by the 7 mm observations. The 
turbulence discussed above accelerates both protons and electrons to a power-law distribution 
with p = 2.2 in the relativistic energy range. The source of electrons for this process is likely 
the hot flow accreting toward the black hole. We therefore assume that the injected electrons 
have a mean Lorentz factor = 10. So the electron distribution has a rising spectrum below 
this energy and a power- law spectrum above it (Park & Petrosian 1995). The distribution 
cuts off when the synchrotron cooling time 

9m^c5 „o / ^ 7 
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equals the acceleration time Tac = 5.0Tesc = 670(4/turb-R^min + 7)(-R/20r5) s. Here, me is the 
electron mass and e is its charge. One therefore gets a high-energy cutoff 

'^'^^ l + 0.6/t,rbi?A:mi„ (tOOg) (20^) ■ ^^^^ 
The corresponding frequency of the synchrotron emission is 

= = [tToIK-^J (looG j (207^) ' ^'^^ 

which should be higher than the observed frequency 43 GHz. We choose I'm > 60 GHz 
(which is 40% above the observed value), giving 

R \ ( 1 \ ( R ^ "^^^ 



S < 126 G < 126 G . (15) 

This constraint is indicated by the vertical hne in Figure 1. 

Because Sagittarius A* has a fiat radio spectrum with a ~ —0.3 {Fy oc z^^"), the optical 
depth of this nonthermal source at 7 mm tq has to be ~ 1 to avoid producing emission in 
excess of the observed flux below (for tq < 1) or above (for tq > 1) 7 mm. We require 
30GHz < < 60GHz, where Tq{vq) = 1 for a source extent ol R — 20r5, i.e. (Pacholczyk 
1970) 

0.08 < f " \ f ^) ( " ( < 0.68 , (16) 
- Vl03cm-3AlO/ VlOOGy \2{)rs)~ 

which excludes the lower-left and upper-right portions of the B — n plane (Fig. 1). The 
luminosity associated with this emission (assuming a spherical source) would be given by 
(Pacholczyk 1970) 



vmL,{vm) = 4.9x10^^ (17) 
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ergs s ^ 



which must be lower than the radio to NIR luminosity of ~ 5 x 10^^ ergs s"^. In combination 
with Equation (11), this constraint excludes the upper portion of the B — n plane (Fig. 1). 
Therefore, to reconcile the radio to IR and 7-ray emissions with the SA model, the plasma 
must be strongly magnetized with /3a > 1- The blank region near the middle of Figure 1 
shows the allowed parameter space, when all these constraints are taken into account. The 
exact properties of the plasma required to explain the 7 mm observations depend on the 
source structure and geometry. As an illustration, we consider a uniform spherical model 
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with radius R — 20 rs (roughly consistent with the VLBA measurement). Calculating the 

accelerated electron spectrum and consequent radio emission yields a best fit model for the 
7 mm data (with B = 70 G, n = 600 cm~^, and /turb-R^min — 0.016) corresponding to the 
black dot near the middle of Figure 1. 

Now, in concordance with the acceleration of electrons to produce Sagittarius A*'s 7 
mm emission, protons are also energized. Figure 2 shows the acceleration and escape time 
scales for electrons {thin) and protons (thick) in plasmas with /3a = 0.15 and 2.0. The time 
scales are given in units of rtj., and we assume /turb-R^min = Svr for this calculation. In 
the general case, Tac oc Tsc oc f^l^ > 1. From these time scales, we see that a hard high 
energy proton spectrum can only be produced in plasmas with /?a > 0.15 and low-energy 
protons have an escape time much longer than their acceleration time. Thus, most of these 
particles are energized efficiently to about E'min ~ 300 MeV, above which the acceleration 
and escape times become comparable and the protons attain a power-law distribution. For 
such a spectrum, we have N{E) = {p — l)n E~p, where we have used the fact that the 
high cutoff energy E'max — ^ireB / kmin — 4.2 x 10^^ eV is much larger than the 'turnover' 
energy £^min- The ensuing energy fiux associated with protons with E > Eg — 0.1 TeV is 
given by 

Fp ' r " (18) 



p-2 T, 



esc 
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With all other parameters being fixed, Fp only depends on n, and is shown on the right axis 
of Figure 1. Note that the illustrative uniform spherical model for the 7 mm observations 
of Sagittarius A* (introduced above) corresponds to a proton energy flux Fp ^ 1.9 x 10^''' 
ergs s~^ The total power associated with the proton flux is ~ 6.0 x 10^'' erg s^^ The 
power injected into the fast wave turbulence can be estimated as ~ (47r/3)i?^VA-B^/turb 
6.0 X lO^^ergs s"^. Most of the dissipated turbulence energy then goes into the protons. And 
for the model parameters O^r < 0.5 for k > /Cmin (note that u ~ E-aAn/lhJneC^ — 60), the 
approximation of the turbulence as parallel propagating waves is therefore appropriate. 



4. The PP Interaction Region 

Once the relativistic protons leave the acceleration site, they diffuse through the medium 
surrounding Sagittarius A*, as described in § 2 above. The transport of cosmic rays through 
chaotic magnetic fields has been considered by several authors, among them Giacalone and 
Jokipii (1999), and Casse, Lemoine, and Pelletier (2001). These particles effectively execute 
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a random walk through the ISM with a mean free path cTsc- There is no evidence of a large 
scale magnetic field on the extension corresponding to the HESS point-spread-function. In 
addition, no structure has been detected within the mini-cavity on angular scales referenced 
to Chandra^s spatial resolution of ~ 1 arcscc. Even more constraining are the numerical 
simulations for the diffuse X-ray emission, which show that shocks resulting from wind- 
wind collisions are typically ten times smaller than the separation between stars. This 
corresponds roughly to a coherence scale ~ 10^^ cm. We therefore infer that /turb — 1 
and /Cmin — 6 X 10"^^ cm~^. Because the temperature of the ambient gas is low (/3p ~ 1), 
the damping of fast mode waves interacting with relativistic protons is unimportant. The 
turbulence is isotropic with g = 3/2 up to which can be estimated with ^cr(^cr) = 1- We 
therefore have kcr — 400A;inin//9p. 

Combined with Equation (5), this therefore yields a scattering time in this region sur- 
rounding Sagittarius A* of r^c 6.4 x 10^(£;/10TeV)^/2(S/0.1 mG)-^/^ s. The HESS source 
may be as large as i? = 3 pc (Aharonian et al. 2004; Aharonian & Neronov 2005b). Thus, 
the total path length of a relativistic proton in the medium surrounding Sagittarius A* is 
/t^t = i?V(rscc) ~ 4.5 X 1022(E/10TeV)-^/2(fi/0.1mG)i/2 cm. On the other hand, the cross 
section (Tpp for pp scattering is weakly dependent on energy, and may be taken as ~ 40 
mb in the energy range of interest (see, e.g., Crocker et al. 2004). Therefore the mean 
free path for a pp scattering event is ~ [(^o)<7pp]~^ ~ 1-3 x 10^^[(no)/200cm~^]~^ cm, so 
hot (no) app ~ 0(1). In other words, most of the escaping relativistic protons scatter with 
background protons about once in the wind-shocked medium surrounding Sagittarius A*. 
Because the corresponding pp scattering time is ~ 10^ yrs, which is longer than the time 
required for the wind- injected ISM to reach equilibrium (Rockefeller et al. 2004), any vari- 
ability in the TeV emission would be associated primarily with density fluctuations on large 
spatial scales, rather than changes in the source. We would expect the source to be stable 
on a time scale of a year or less.^ 

The particle cascade generated by scatterings between the relativistic protons diffusing 
away from the acceleration site and protons in the ambient medium produces additional 
spectral components due to neutral pion decays and emission by secondary leptons. Neutral 
pion decays produce a 7-ray spectrum that dominates above ~ 100 MeV. At lower energies, 
most of the radiative emission is due to the decay products of charged pions. These leptons 
radiate via synchrotron, bremsstrahlung, and inverse Compton scattering with the strong 
UV and IR fields bathing the inner few parsecs of the Galaxy (Telesco, Davidson & Werner 



^Equation (6) shows that the corresponding acceleration time in this region is t^c — 7.5 x 

10^Tsc((no)/200cm^'^)(iJ/0.1 rriG)^^, which is longer than any other relevant time scales, so the acceler- 
ation of protons by the turbulence outside of Sagittarius A* itself may be ignored. 
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1996; Davidson et al. 1992; Becklin, Gatley & Werner 1982). A detailed description of 
how the total emission from these various processes can be calculated appears in Fatuzzo 

6 Melia (2003). For a proton energy flux Fp ^ 1.6 x 10'^'' ergs s~^ (above 300 MeV with 
p = 2.25), the broadband spectrum resulting from these processes is shown in Figure 3. 
For this calculation, we assume that the protons are scattered only once and the magnetic 
field strength in the pp interaction zone was set at 1.0 mG. The currently available data 
(or upper limits) are also shown here for comparison with our model. We also include the 
GLAST single energy sensitivity for a one-year all-sky survey, which shows that the HESS 
~ 100 MeV-100 GeV counterpart may also be detectable. Note that since coohng time of 
TeV electrons and positrons in the pp interaction region is considerably shorter than the pp 
scattering time (~ 4 x 10^^ s), their distributions reach a cooling broken-power spectrum in 
the steady-state. 

The spectrum fits the HESS data and is consistent with observations in other lower 
energy bands. The fact that the required proton power is more than 3 times below that 
predicted from the fitting of the 7 mm data (~ 6.0 x 10^'' ergs s~^) indicates that most 
of the protons diffuse beyond the HESS point- spread-function, which may explain the re- 
cently observed diffusion gamma-ray emission from the Galactic Center (Aharonian et al. 
2006), constituting one of the strongest arguments in favor of the SA model for the origin of 
Sagittarius A*'s spectrum. 

5. Conclusion and Discussion 

The combination of a reliable intrinsic source size measurement of Sagittarius A* at 

7 mm with the detection of TeV emission from its vicinity has important implications for 
the acceleration of relativistic particles near the black hole. We have shown that SA of 
electrons and protons in a magnetic field dominated plasma can account for both the 7 mm 
observations via in situ electron synchrotron process and the HESS observations through 
the diffusion of accelerated protons into the surrounding ISM. The fact that a medium with 
/3a > 1 is required is intriguing in fight of recent MHD simulations of black hole accretion 
(Hawley & Balbus 2002; Hirosc ct al. 2004). These numerical calculations show that there 
are always magnetic-field-dominatcd coronas above the disk and funnels along its symmetry 
axis. These regions may be the electron and proton acceleration sites. Speculating further, 
our results may also be in line with current concepts of a Poynting flux dominated outflow, 
driven either by an accretion torus or by a spinning black hole (Camenzind 2004; Semenov, 
Dyadechkin, & Punsly 2004; Nakamura & Meier 2004). 

Although we suggested that most of the protons producing the TeV emission scatter 
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at least once within a region corresponding to the HESS point-spread-function, given the 
uncertainties in the observed properties of the ISM, it is still possible that a significant 
fraction of very high energy protons diffuse outward to even larger radii. Then the proton 
spectral index in the acceleration region would need to be smaller than 2.2 in order to explain 
the HESS photon spectrum since the proton scattering time increases with energy. However, 
the difference between the proton and photon spectral indexes, which depends on the energy- 
dependent diffusion of the former, is expected to be small (Aharonian & Neronov 2005b), 
and our main conclusions still hold. In this paper we studied the SA of particles interacting 
with parallel propagating waves. This may not be true at the wave injection scale where the 
turbulence could be isotropic. A comprehensive investigation of this aspect will be presented 
in a separate paper. 

We have shown above that particles energized in a magnetized plasma apparently pro- 
duce most of Sagittarius A*'s emission. This may very well be the case for the accretion 
flow in many other compact objects. A thorough investigation of the SA model for parti- 
cles in a magnetic field dominated outflow will also play an important role in bridging the 
above mentioned theoretical developments with observations. Eventually, studying SA in 
such compact environments may also reveal the nature of large scale relativistic jets, where 
the coupling of plasma turbulence to the energetic particles could be very strong. 

The radiative consequences of relativistic protons in active galactic nuclei and their 
implications on the formation of large scale jets have been studied extensively during the 
past few decades (Bcgclman, Rudak, & Sikora 1990; Atoyan & Dcrmcr 2003). With the 
theory of SA developed here, properties of the turbulent plasma responsible for the proton 
acceleration can be inferred. The constraints to the particle acceleration processes would be 
stricter should the acceleration of electrons by the same turbulence be also considered, as is 
the approach we took here in modeling the broadband spectrum of Sagittarius A*. 

We have used the spatially integrated results of SA as a first order approximation in 
the previous investigations. However, to model spatially resolved sources, one must take 
into account the source structure. Monte Carlo simulations of SA may be a good approach 
in such cases and this will hkely be the next step in the development of our model. For 
a complete treatment of SA, the properties of the waves must also be addressed. MHD 
simulations provide a complementary vehicle for this problem and when combined with 
the SA theory, can yield a global picture of energy dissipation channeled via turbulence. 
Investigations of this kind will be fruitful in providing a better understanding of many 
nonthermal astrophysical sources. With little modification, the model can also be applied to 
other astrophysical systems where the plasma is coUisionless and strong turbulence presents. 
The observed emission spectra can be used to determine the properties of the plasma and 
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to reveal the underlying physics processes. 
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log,o(B/Gauss) 

Fig. 1. — Constraints on the density and magnetic field of the acceleration site, where 
electrons produce the observed 7 mm emission while protons escape toward large radii and 
produce the HESS 7-ray signal via collisions with the ambient protons. Here we have assumed 
that 7in = 10 and R = 20 r^. The dashed line corresponds to (3a = 1- At 7 mm the source 
is optically thin (thick) in the lower-left (upper-right) portion of the figure. The constraint 
S < 120 G comes from the requirement that the synchrotron emission cuts off above 60 
GHz. The upper area is excluded because the total synchrotron power from the acceleration 
site can not exceed the radio to IR luminosity of Sagittarius A*. Depending on the details 
of the source geometry, the 7 mm observations may be explained with B and n in the blank 
(central) region of this figure. For example, the black dot corresponds to an illustrative 
uniform sphere model. The right axis indicates the power in escaping protons with energy 
E >Q.l TeV. 
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Fig. 2. — The acceleration (solid) and escape (dashed) time scales for electrons (thin) and 
protons (thick), in plasmas with (3 a = 0.15 and 2.0 (indicated in the figure). The time 
scales are shown in units of the hght transit time rtr = R/c. The turbulence is assumed 
to be in energy equipartition with the magnetic field and q = 2. In a weakly magnetized 
plasma (/3a < 1-0), whistler waves can accelerate low-energy electrons to a Lorentz factor 
13 Atrip /mei causing a sharp rise in the electron acceleration time with energy beyond 1 MeV. 
The acceleration of nonrelativistic electrons is dominated by high frequency electromagnetic 
waves. In a strongly magnetized plasma {(3a > 1-0), the whistler wave branch vanishes, the 
acceleration time of relativistic electrons is almost energy-independent. 
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Fig. 3. — Broadband emission from the pp interaction zone due to (short dashed curve) 
neutral pion decays, and (long dashed curve) bremsstrahlung, (short dash-dot curve) UV in- 
verse Compton scattering, (long dash-dot curve) IR inverse Compton scattering, and (dotted 
curve) synchrotron emission by the secondary leptons. The solid curve represents the total 
broadband emission. For comparison, we also show the currently available data (or upper 
limits), and the expected GLAST single energy sensitivity for a one-year all-sky survey The 
Integral observations arc from Bclangcr ct al. (2005), the EGRET source 3EG J1746-2852 
is from Mayer-Hasselwander et al. (1998), and the HESS observations are reported in Aha- 
ronian et al. (2004). The Chandra upper limit may be found in Rockefeller et al. (2004), 
and the VLA datum is reproduced in Fatuzzo and Melia (2003). 



